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Abstract 

In this paper we consider the two atoms Tavis-Cummings model and give an 
exphcit form to the solution of this model which will play a central role in quantum 
computation based on atoms of laser-cooled and trapped linearly in a cavity. 

We also present a problem of three atoms Tavis-Cummings model which is re- 
lated to the construction of controlled-controlled NOT operation (gate) in quantum 
computation. 
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The purpose of this paper is to give an exphcit form to the solution of Tavis-Cummings 
model (P) with one and two atoms. This model is a very important one in Quantum 
Optics and has been studied widely , see |21, |S] or |lj as general textbooks in quantum 
optics. See also recent papers 0, [Hj and their references. 

We are studying a quantum computation and therefore want to study the model from 
this point of view, namely the quantum computation based on atoms of laser-cooled and 
trapped linearly in a cavity. We must in this model construct a controlled NOT gate or 
other controlled unitary gates to perform a quantum computation, see [7] as a general 
introduction to this subject. 

For that purpose we need the explicit form of solution of the models with one, two 
and three atoms. As for the model of one atom it is more or less well-known, and as for 
the case of two or three atoms it has not been given as far as we know. 

In this paper we give it for the case of two atoms, while we could not give it for the 
three atoms case, so we present it as a challenging problem. Anyway, let us start. 

The Tavis-Cummings model (with n-atoms) that we will treat in this paper can be 
written as follows (we set h = 1 for simplicity). 



^ n n 

H = ujlL0a^a + — Y^ (xf ^ ® 1 + ^ ^ (g) a + aj"^ ® 



2 i=i i=i 
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where uj is the frequency of radiation field, A the energy difference of two level atoms, 
a and are annihilation and creation operators of the field, and g a coupling constant. 



and L = 2". Here al , al and al are given as 



al'> = I2 ® • • ■ ® I2 ® ® I2 ® ■ ■ ■ ® I2 {i- position) e M{L, C) (2) 



where s is +, — and 3 respectively and 

l^ / 

0"+ = I h 0"- = 1,0-3 






Here let us rewrite the hamiltonian If we set 



1 
VO -1 



(1 

U 1 



(3) 
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then can be written as 

H = ujIl ® a'^a + AS3 ^ 1 + g (^S+ ® a + S_ a'^^ =Ho + V, (5) 

which is very clear. We note that {S^, S^, S3} satisfy the sM(2)-relation 

[53, 5+] = [53, S^] = -5_, [5+, 5_] = 253. (6) 

However, the representation p defined by 

p(o-+) = 5+, p{a^) = S^, p{as/2) = S3 

is a reducible representation of su{2). 

We would like to solve the Schrodinger equation 

tjU = HU= {Ho + V) U, (7) 

where [/ is a unitary operator. We can solve this equation by using the method of 
constant variation. Let us make a brief review. The equation i-^U = HqU is easily 
solved to be 

U{t) = (e-^*^^^ ® e"'*"^) f/o 

where N = a^a is a number operator and Uq a constant unitary. By changing Uq i — > 
Uo{t) and substituting into ((Zj) we have the equation 

i^Uo = g {e'*(^-")5+ ® a + e-'*(^-")5_ ® a^} Uo (8) 

after some algebras. Here let us assume the resonance condition 

A = (9) 

which makes the situation simpler. Under this condition the solution of (jHl) becomes 

[/o(t) =e-**^(^+®"+^-®''^), (10) 

so that the full solution of ((Tj) is given by 

U{t) = (e-**'"'^^ ® g~itc.7V^ ^-itg{S+^a+S.<^at) ^ ^^^^ 

where we have dropped the constant unitary operator for simplicity. Therefore we have 

only to calculate the term explicitly, which is however a very hard task ^. In the 
^the situation is very similar to that of 
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following we set 



A = S+ ® a + (g) 



'12) 



for simplicity. We can determine e ^'^"^ for n = 1 (one atom case) and n = 2 (two atoms 
case) completely. Let us show. 




One Atom Case In this case A in ()12p is written as 



A 



Since 

aa'' 

a'^a 

with the number operator and 



A' 



N + l 
\ N 



D 



A^3 = D\ A^-'+i = D^A, for j > 0, 



(13) 



(14) 



so we have 



'itgA 



t^} ^ 27 ! ^ 27 + 1 ! 



j=0 

oo 

V(-l)^ 



j=0 



j=0 



(2j + l)! 



cos 



.sin(igV]v+T) 



V 



.sinftg^Jv) X 
—i — -a^ 



(15) 



" Vn+1 ^ 

We obtained the explicit form of solution. However, this form is more or less well-known, 
see for example Chapter 14. Because this solution is very convenient, there are many 
applications, see the textbook 0. 

We note that (|TH|) can be decomposed as 



/ 



cos 



.sm(tg^/N+T] 



.SmitgVN) . 



VN+1 

COS (tgVN^ 



cos 



-la 



V7V+T 



.Sm.{tg^/N+l) 

cos (tgVN^ 



V 

/ 

V 



-la 



1 

^tdJl(tgy/N+l) 



( 



V 



cos 







1 



—I 







tan(igv^) J. 

(X -L 



1 



COS [tg^jN + 1 




COS(t9Vlv) / 

^ / 




COS(igVJV) / 



1 



.tan(tgV]v+T) 



\/]v+T 
1 



(16) 



We leave the check to the readers. 

This is a Gauss decomposition of unitary operator. This may be used to construct 
a theory of "quantum" representation of a non-commutative group, which is now under 
consideration. 

Two Atoms Case In this case A in ()12|) is written as 



A 



a a ^ 

at a 

at a 

at at 



(17) 



We would like to look for the explicit form of solution like (jlSp . so we must find a relation 
like (fT^ . It is not difficult to see 



A^ 



^ 2(A^+ 1) 2a2 ^ 

2A^ + 1 2A^ + 1 

2A^ + 1 2A^ + 1 

2(at)2 ^ IN j 



(18) 



A' 



(2iV + 3)a (2A^ + 3)a 
(2A^+l)at 
(2A^+l)at 

(2A^-l)at (2Ar-l)at 





(2iV + l)a 
(2iV + l)a 




(19) 
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From ()19p we find a clear relation 



A' 



2(2A^ + 3) 



2{2N +1) 



2{2N + 1) 



2(2iV- 1) ) 



^ ^ a a ^ 



at a 
at a 
at at 



DA. 



This is our key observation. From this it is easy to see 



(20) 



SO that we have 



-itgA 



1 + E 



oo 



^^d]^} m ^ ^VlJ^y ^) (2, + !)! 



1 + {-1 + COS {tgVD^ }A^- ^-^sin {tg\/l5^ A 



(21) 



or more explicitly 



-itgA 



^ ail 0,12 0-13 (3.14 ^ 



a21 a22 0,23 (^2A 
C^Sl '^32 a34 

y a4i a42 043 a44 j 



where 



(22) 



N + 2 + {N + l)cos (t^y^2(2A^ + 3)^ 

flu = „ , „ , fll2 = fll3 



sin (t^y^2(2iV + 3)) 



2iV + 3 



2(2A^ + 3) 



ai4 



-1 + cos (tg^2{2N + 3 
2ivT3 



a , a2i — a3i 



sin t^J2(2iV+ 1) 



2(2iV + 1) 
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^22 — 0-33 



1 + COS (t^^2(2A^+ 1)) 



^44 



-1 + COS (tg^2{2N -I) 

2N -I ' 

- 1 + A^cos (tg^2{2N -I' 
2N -I 



4^2 



(^23 — '^32 



042 — Cf43 



-1 + COS (tg^2{2N + 1 



sin [tgJ2{2N -I) 



2(2iV- 1) 



This is our main result in this paper. The exphcit form has not been known in the 
hterature as far as we know. 

Since the Tavis-Cummings model has a kind of universal characteristic and the ex- 
plicit form of solution was given, there must be many applications to Quantum Optics, 
Mathematical Physics and etc. In the forthcoming paper J2] we will apply this to the 
construction of controlled-unitary operations (gates) in quantum computation (see for 
example [Zj) based on atoms of laser-cooled and trapped linearly in a cavity. 

Three Atoms Case In this case A in (fT^ is written as 

'^OaaOaOOO^ 
atOOaOaOO 
atOOaOOaO 
at at a 
atOOOOaaO 
at at a 
at at a 
at at at 



A 




(23) 



where A is A in (fTTj) . 

We would like to look for the explicit form of solution like (|T5|) or ()22|) . However, 
we could not find a relation like (fT^ or ()2()|1 for ()23|1 in spite of much effort (we have 
calculated A^, ■ ■ -, A^). We encourage the readers to tackle this problem. 

We note that the solution is deeply related to the construction of controlled-controUed 
NOT operation (gate) in quantum computation, so the explicit form of it is needed. 



We conclude this paper by making a comment on our target. The Tavis-Cummings 
model is based on two energy levels of atoms. However, an atom has in general infinitely 
many energy levels, so it is natural to use this possibility. We are studying a quantum 
computation based on multi-level systems of atoms (a qudit theory), lEj^PB]- Therefore 
we would like to extend the Tavis-Cummings model based on two-levels to a model based 
on multi-levels. This is a very challenging task ! 



Appendix 



In this appendix we show an another approach to obtain the result in the two atoms case 
which may be useful in the three atoms one. Our method is to reduce the 4 x 4-matrix 
A in ()17|1 to a 3 X 3-matrix B in to make our calculation easier. For that aim we 
prepare the following two matrices 



/ 



T 



1 



1 

V2 

1 
V2 



1 

1 
V2 



\ 



s 



( 



1 



Then it is easy to see 



s (tat^) = (ST) A {sry 













V2at 







V2a^ 










B 



where 



B 





v^at 




^/2a 


V2a^ 








J+ ® a + J_ 



(24) 
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{ J+, J_} are just generators of (spin one) irreducible representation of Therefore to 
calculate e~**^^ we have only to do e~**^'^. The method is almost similar. Namely, 



= 2 



N + 1 

2A^ + 1 

t2 



a' 



N 



I 



2(2iV + 3) 



2(2A^ + 1) 



2(2iV-l) ] 



B = DB, 



so we obtain 



Q-itgB _ 1 + ^ |_i ^ cog (^tgVD)}B^ - ^-^sin {igVo' 



B 



^ bn hi2 bis ^ 

^21 ^22 ^23 
\ hi &32 ^33 J 



(25) 



where 



(^11 

b22 

bsi 

&33 



+ 2 + (A^ + l)cos (tg^2{2N + 3)) sin (t^^2(2A^ + 3) 

! bi2 = —i— 



2N + 3 



1 + cos (tg J 2{2N + 3) 



2N + 3 



a , b 



21 



V2iVT3 

sin (tgJ2{2N+l)) 

^ — , -a^ , 

^2N+l 



cos [tg J 2{2N + I)], h 



J23 



-1 + cos {tg^2{2N-l)) ^ 

2N-1 ' 



sin (tg^2{2N+l) 
V2N + 1 

sin (tg^2{2N - I' 



V2iV- 1 



iV - 1 + iVcos (t^^2(2iV - 1)) 



2A^- 1 



We leave the remainder to the readers. 
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Note added in the text 

After submitting this paper Pablo P. Munhoz kindky informed us of the paper 
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M.S. Kim, J. Lee, D.Ahn, and P.L. Knight : Entanglement induced by a single-mode heat 
environment, Phys. Rev. A 65, 040101 (2002), |quant^ /0109052| 

and suggested seeing the formula (4). We have checked the agreement between the formula 
(4) and our result, so our result is not new. However, how to derive the formula is not 
written in the paper and our method seems much simpler (or easy to understand), so our 
paper is still valuable. 

We wish to thank Pablo P. Munhoz for his helpful comment. 
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